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Charge and orbital ordering in the low-temperature monoclinic structure of magnetite (Fesd) 
is investigated using LSDA+t/. While the difference between t2 g minority occupancies of Fe^ + and 
Fe 3 ^ 1 " cations is large and gives direct evidence for charge ordering, the screening is so effective 
that the total 3d charge disproportion is rather small. The charge order has a pronounced [001] 
modulation, which is incompatible with the Anderson criterion. The orbital order agrees with the 
Kugel-Khomskii theory. 

PACS numbers: 75.10 Jm, 75.30.Gw, 75.70.Ak 



The magnetic properties of magnetite, the famous 
lodestone, has fascinated mankind for several thousand 
years already 0]. Even today, in view of the possible 
technological importance of this material for spintron- 
ics 0, and because of the still not well understood low- 
temperature properties, magnetite remains at the focus 
of active research. 

Magnetite is a mixed-valent system and is the par- 
ent compound for magnetic materials such as maghemite 
(Fe2C>3) and spinel ferrites. At room temperature it crys- 
tallizes in the inverted cubic spinel structure Fd3m with 
tetrahedral A-sites occupied by Fe 3+ cations, whereas oc- 
tahedral _B-sites are occupied by an equal number of ran- 
domly distributed Fe 2+ and Fe 3+ cations. Magnetite is 
fcrrimagnetic with an anomalous high critical temper- 
ature Tc ~ 860 K, the ^4-site magnetic moments being 
aligned antiparallcl to the i?-site moments. At room tem- 
perature Fe3C>4 is a poor metal with an electronic con- 
ductivity of 4 mfi cm. Upon further cooling a first-order 
metal-insulator (Verwey-) transition occurs at Ty ~ 120 
K where conductivity abruptly decreases by two orders of 
magnitude. According to Verwey this transition is caused 
by the ordering of Fe 2+ cations on the B sublattice, with 
a simple charge arrangement of (001 ) c planes (indexed on 
the cubic cell) alternately occupied by 2+ and 3+ Fcb 
cations (Verwey charge ordering model) |3j, |4| . This par- 
ticular charge order (CO) obeys the so-called Anderson 
criterion [j| for minimal electrostatic repulsion leading to 
a short range CO pattern, namely tetrahedra of B-sites 
with an equal number of 2+ and 3+ cations. Since then a 
wide range of other CO models has been proposed which, 
however, all make use of the Anderson criterion 



Later experiments showed that the Verwey transition 
is accompanied by a structural distortion from cubic to 
the monoclinic structure which has not been fully re- 
solved so far 0, Q • Although the experiments rule out 
the Verwey CO below Ty all attempts to construct a re- 
fined CO model from neutron-scattering data set failed 
because neutron scattering is more sensitive to atomic 



displacements than to charge ordering. In the absence of 
a definitive, experimentally determined structure many 
theoretical models for the low-temperature (LT) phase of 
magnetite |9J[have been proposed. They include purely 
electronic 0, Q] and electron-phonon 0, 0] models 
for CO, as well as a bond dimerized ground state with- 
out charge separation v 14jj. In particular all previous 
LDA(+[/) calculations were performed for undistorted 
cubic unit cell. In spite of the fact that the amplitudes 
of these distortions are quite small this approximation 
for the LT unit cell in LSDA+t/ calculations results in- 
evitably in the Verwey CO. This problem is overcome in 
our work using recently refined crystal structure [TH fl6| 
in which the ground state with more complicated CO 
was found. This also confirms Szotek et al. [TtI ] con- 
clusion that the Verwey CO is not the ground state for 
magnetite, which is based on self-interaction corrected 
local spin density calculations for the refined LT crystal 
structure. 

In this paper we report LSDA+C/ calculations [Til llflj 
in the tight-binding linear muffin-tin orbital (TBLMTO) 
calculation scheme for Fe304 in the P2/c structure 
[Tfl flflj]. Motivated by our results, we propose an order 
parameter, defined as the difference between t2 g minority 
spin occupancies of Fe 2 ^ and Fe 3 ^ cations. This order 
parameter is found to be quite large, although the total 
3d charge difference between these cations, is small. It 
seems certain that magnetite is long-range ordered be- 
low Ty, in contrast to the intermediate valence regime 
proposed by Garcia et al. [2(1 EH E^ . 

Recently, the LT superstructure of magnetite was re- 
fined by Wright et al. [lil flfij ]. The space group was 
confirmed to be monoclinic Cc, but the structure refine- 
ment was only possible in P2/c group. They proposed 
the a c /\/2 x a c /y2 x 2a c subcell (a c is the cell parame- 
ter of the undistorted cubic unit cell) with P2/c space 
group symmetry. The additional Pmca orthorhombic 
symmetry constrains were also applied. The refined cell 
parameters were a = 5.94437(1) A, b = 5.92471(2) A, 
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c = 16.77512(4) A, and (3 = 90.236(1)°. The structural 
evidence for the CO below the transition in the refined 
crystal structure, which is based on estimations of the 
mean 5-site-to-oxygen distances or the bond valence sum 
(BVS) analyses, was also declared 0,0. However, this 
refined structure analysis has recently been found to be 
controversial |2l|. The lack of atomic long-range order 
and, as a result, the intermediate valence regime below 
the Verwey transition were proposed. Indeed, a differ- 
ence between averaged Fe-0 distances of compressed and 
expanded FeOg octahedra corresponds to the maximum 
limit of charge disproportion (0.2e), which has the same 
order as the total sensitivity (including experimental er- 
rors) of BVS method. This contradiction as well as the 
ambiguity of proposed CO schemes (two different CO 
classes were proposed: class-I in the P2/c unit cell and 
class-II in the full Cc superstructure), is resolved in our 
electronic structure study. 

We perform LSD A and LSDA+[/ calculations for 
Fes04 in the P2/c structure. For simplicity we neglect 
the small spin-orbit coupling (in previous calculations 
for cubic Fe304 spin-orbital splitting of the 3d-band was 
found to be two orders of magnitude smaller than the 
crystal field splitting HiJ). The LSD A calculations give 
only a half-metallic ferrimagnetic solution without CO. 
Partially filled bands at the Fermi level originate from 
the minority spin t2 g orbitals of Fes cations. The tetra- 
hedral Fe^ cations do not participate in the formation 
of bands near the Fermi level, since their minority and 
majority spin 3d states are completely occupied and com- 
pletely empty, respectively. Thus, the LSDA results qual- 
itatively agree with previous band-structure calculations 
for the cubic phase |rjj,|24j,|25j, 26]. Apparently, only crys- 
tal structure distortion from cubic to monoclinic phase 
is not sufficient to explain metal-insulator transition and 
charge ordering in magnetite. The electron-electron cor- 
relations, mainly in the 3d shell of Fe cations, play sig- 
nificant role. 

To proceed further we take into account the strong 
electronic correlations in Fe 3d shell using LSDA+i7 
method. In contrast to the LSDA, already with U and 
J obtained from constrained calculations (4.5 eV and 1 
eV, respectively) a charge ordered insulator with energy 
gap of 0.03 eV was obtained. On the other hand, the 
calculation of U depends on theoretical approximations 
and, as a rule, the accuracy does not exceed 10-20%. 
A reasonably good agreement of the calculated gap of 
0.18 eV with the experimental value of 0.14 eV |2jJ was 
found using the U value of 5 eV. As shown in Fig. ^ the 
gap opens between occupied and unoccupied t^gi states 
of Fe^,, Fe^ 4 ^ and Fe^t^, Fe^"^ respectively. The 
remaining unoccupied Fes states are pushed by strong 
Coulomb repulsion to the energies above 2 eV. 

The obtained CO is described by a dominant [001] c 
charge density wave with a minor [00^] c modulation in 
the phase of CO and coincides with the class-I CO, which 
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FIG. 1: Partial DOS obtained from the LSDA+C7 calcula- 
tions with (7=5 eV and J=l eV for the low-temperature P2/c 
phase of Fe304. The top of the valence band is shown by dot- 
ted lines. 



was proposed earlier by Wright et al. Thus, the 

LSDA+C/ calculations confirm violation of the Anderson 
criterion for Fe304 in the LT phase. In order to check the 
stability of the obtained CO solution we performed addi- 
tional self-consistent LSDA+t/ calculations both for the 
P2 1 c structure and Cc supercell. In the first case the Ver- 
wey CO was used as the starting CO model, while in the 
second one we started from class-II CO, shown in Fig. 2 
of Ref. 16j However, it was found that these CO models 
are unstable and the obtained self-consistent solutions 
coincide with the stable one found previously, i.e. for 
certain value of U the obtained CO does not depend on 
initial charge arrangement. Obviously, we did not check 
all possible CO scenarios within P2/c unit cell or Cc 
supercell, but our results consistently indicate that the 
obtained class-I CO solution is the ground state of Fea04 
in the LT phase. It is important to note that LSDA+C7 
calculations performed for an undistorted P2 / c supercell 
of Fd3m structure result in an insulating CO solution 
which is compatible with the Verwey CO model. 

An analysis of occupation matrices of 3<i-Fe b minority 
spin states confirms substantial charge disproportiona- 
tion. As shown in Tabled one of the states of Fe^^ 

and Fe^t^ cations is almost completely filled with the 
occupation numbers n w 0.8. On the other hand, the 
remained two t2gi orbitals of the Fe 2+ cations have sig- 
nificantly smaller population of about 0.15. The occu- 
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TABLE I: Total and /-projected charges, magnetic moments, 
and the occupation of the most populated t^g minority orbital 
calculated for inequivalent Fes ions in the low-temperature 
P2/c phase of Fe 3 4 [H. 



Fes ion 


q 


Is 


q P 


Id 


M (>b) 


t2 g [ orbital 


n 


Fe B(1) 


6.04 


0.17 


0.19 


5.69 


3.50 


dxz -F dy z 


0.76 


Fe B(2) 


5.73 


0.19 


0.21 


5.44 


3.94 




0.27 


Fes (3) 


5.91 


0.19 


0.21 


5.51 


3.81 




0.24 


Fe S ( 4 ) 


6.03 


0.16 


0.18 


5.69 


3.48 


(i X 2 y 2 


0.80 



pation numbers of t 2g i orbitals for Fe^ 2 ^ and Fe^ 3 ^ do 
not exceed 0.3, which gives the value of about 0.5 for 
the largest difference between Fe^ + and Fe^ + t^g minor- 
ity spin populations. The corresponding total 3d charges 
difference (0.23) and disproportionation of the total elec- 
tron charges inside the atomic spheres of Fe^ + and Fe 3 B + 
cations (0.24) is in a reasonably good agreement with the 
value of 0.2 estimated from BVS analysis of monoclinic 
structure. This shows that the change of the t2 3 i occu- 
pations caused by the charge ordering is very effectively 
screened by the rearrangement of the other Fe electrons. 

Significant contribution to the charge screening is pro- 
vided by Fes e g states. Although the bands originating 
from these states are located well above the energy gap, 
the e g minority orbitals form relatively strong tr-bonds 
with 2p-states of the oxygen octahedron and, as a result, 
give appreciable contribution to the occupied part of the 
valence band. The above mentioned screening reduces 
the energy loss due to the development of CO incompati- 
ble with the Anderson criterion in the LT phase of Fe3C>4. 
Hence, due to strong screening effects, the order parame- 
ter, introduced as the total 3d charge difference between 
2+ and 3+ Fes cations, is ill-defined. That explains why 
the BVS analysis does not give convincing proof of CO 
existence. Apparently, the well-defined order parameter 
is the difference of t 2g i occupancies for Fe^ + and Fe^ + 
ions, which amounts to 50% of ideal ionic CO model and 
clearly pronounces the existence of CO below the Verwey 
transition. 

The self-consistent solution obtained from the 
LSDA+U calculations is not only charge but also or- 
bitally ordered. The occupied minority spin tig state of 
Fe^tjj and Fe^ 4 ^ cations is predominantly of d xz ± d yz 
and d x 2_ y 2 character, respectively. This is illustrated by 
Fig. which shows the angular distribution of the mi- 
nority spin 3d Fes electron density. Note, however, that 
the P2/c frame is rotated by an angle ~ 7r/4 with re- 
spect to the cubic one and the angular dependence of t 2g 
states is given by d xy ± d yz and d x 2_ y 2 combination of 
cubic harmonics. The obtained relative orientation of oc- 
cupied Fes tig minority orbitals corresponds to the anti- 
ferro-orbital order. Since all Fes cations are ferromag- 
netically coupled, the obtained orbital order is consistent 
with the anti-ferro-orbital ferromagnetic state, which is 




FIG. 2: The LSDA+[/ angular distribution of the 
minority spin 3d electron density of Fes cations for 
the low-temperature P2/c phase of Fes04. The an- 
gular distribution is calculated according to p(8, (f>) = 
2~2 m ,m' n m, m ' Y^(8 , <f>) Y m t(0,</>), where n m>m / is the occu- 
pation matrix of d minority states for Fes atoms. Y m (9,4>) 
are corresponding spherical harmonics. Oxygen atoms are 
shown by small spheres. The size of orbital corresponds to its 
occupancy. 



the ground state of the degenerate Hubbard model ac- 
cording to the Kugel-Khomskii theory [29|. This orbital 
order leads to the corresponding distortions of FeOa octa- 
hedra. An analysis of interatomic distances in the mono- 
clinic structure ( Table ITTf) shows that the average Fes(i)- 
O distance (2.087 A) in the plane perpendicular to one 
of the diagonals of the distorted FeB(i)Oe octahedron is 
considerably larger than average distances in the other 
two planes (2.063 and 2.067 A). It turns out that the 
occupied FeB(i) t 2g minority spin orbital is the one ori- 
ented in the plane with the largest average Fq-b(\) O dis- 
tance. The same is also true for Fe^(4) ion but in this 
case the variation of the average Fe^^j-O distances is 
smaller (2.074 A vs. 2 x 2.067 A) and, as a consequence, 
the out-of-plane rotation of the occupied t 2g [ orbital is 
stronger. 

As was shown earlier, the Verwey CO model possesses 
the minimum electrostatic repulsion energy among all 
possible CO models But due to the existence of 

two perpendicular families of B-site chains (for instance 
[110] c and [110] c ) correspondingly occupied by 2+ and 
3+ Fes cations the lattice "feels" significant stresses and 
tends to expand in one ([110] c ) and to compress in the an- 
other ([110] c ) direction. Therefore, the Verwey CO gives 
a significant "elastic" energy contribution into the to- 
tal energy and in spite of the lowest electrostatic energy 
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TABLE II: The averaged Fes-0 distances in the plane of tig 
orbitals for P2/c structure of Fea04 (2^1 



Fes atom 


orbital 


d orb . (A) 


d av . (A) 


Fee (ia) 


r\ 4-/7 

!*XZ \ U>yz 




9 079 


d,XZ dyz 


2.087 






d x 2 —y 2 


2.063 




Fe B (ib) 


dxZ ~\~ dyz 


2.087 


2.072 


dxZ ' dyz 


2.067 






dx 2 — y 2 


2.063 




Fe B ( 4 ) 


dxZ i dyz 


2.067 


2.069 


d x 2 —y 2 


2.074 





becomes less favorable than other arrangements. The 
competition of these two (elastic and electrostatic) con- 
tributions in the total energy appears to be responsible 
for the charge order, which is realized in the experimen- 
tally observed low-temperature monoclinic structure. In 
this CO scheme the alternating (001) c planes occupied 
by 2+ ("occupied" plane) and by 3+ ("empty" plane) 
Fes cations are separated by the "partially" occupied 
plane. This (001) c planes order makes the difference be- 
tween [110] c and [110] c directions less pronounced and 
significantly reduces the lattice stress and, as a result, 
reduces the elastic energy contribution in the total en- 
ergy. We propose that this scenario is the primary cause 
for development of the class-I CO found below the Ver- 
wey transition. 

In summary, in the present LSDA+J7 study of the LT 
P2/c phase of Fe304 we found a charge and orbitally 
ordered insulator with an energy gap of 0.18 eV. This 
is in a good agreement with the experimental value of 
0.14 eV [27j. The obtained charge ordered ground state 
is described by a dominant [001] c charge density wave 
with a minor [00^] c modulation. The CO coincides with 
the earlier proposed class-I CO (l ajlq and confirms a 
violation of the Anderson criterion [5j ■ While the screen- 
ing of the charge disproportion is so effective that the 
total 3d charge disproportion is rather small (0.23), the 
charge order is well pronounced with an order parameter 
defined as a difference of t2 9 i occupancies of 2+ and 3+ 
Fes cations (0.5). This agrees well with the result of BVS 
analysis for monoclinic structure (0.2). The orbital order 
is in agreement with the Kugel-Khomskii theory [2^| and 
corresponds to the local distortions of oxygen octahedra 
surrounding Fes-sites. 
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